Binding of rj c meson with light nucleus. 
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Abstract 

Binding energy of r\ c — 4 He system has been calculated on the basis of Yukawa 
type r] c — N potential, suggested in [1]. Final Rank Approximation (FRA) has been 
used for microscopical treatment of the initial 5-body problem. The results are 
compared with the results of the folding model and variational calculations. 

As pointed out in [1] there is an interesting peculiarity in the rj c — N interaction. It 
turned out that the main contribution to the meson-nucleon coupling is due to multi- 
gluon exchanges, a property modelled by a simple Yukawa potential. It means that via 
the study of nonrelativistic systems like rj c plus few nucleons one can get information 
about the role of the gluonic degrees of freedom. 

Previous considerations [2] of rj c - meson interaction with light nuclei (A = 3,4) were 
based on the so called folding model. In this paper, starting from four- and five- body 
meson-nuclear Hamiltonians, we will develop a more elaborate scheme which contains the 
folding model particular case. 

Let us write the total Hamiltonian for the system of one meson and A nucleons in the 
form 

A 

H = h + J2 Vi + H A = ho + V + H A . (1) 

t=i 

Here ho denotes the kinetic energy of the meson relative to the center of mass of the 
nuclear subsystem, Vi is the interaction of the meson with the i-th nucleon, and H A the 
Hamiltonian of the nuclear subsystem. 

In what follows we are interested only in calculating the binding energy of meson 
with 4 He. Having in mind, that in such a situation we are comparatively far away from 
the nearest break-up thresholds for nuclear subsystem, we can approximate in eq.(l) the 
nuclear subsystem Hamiltonian H A by a finite rank operator, for instance, by the rank-one 
term 

H A ^ H A = e \xo)(Xol (2) 
where e an d Xo denote the ground state energy and the corresponding wave function of 
3 He or 4 ife respectively. This Finite Rank Approximation (FRA) is more general, than 
the Resonating Group Method (RGM), as will be seen below from the expression for the 
eigenf unctions of the Hamiltonian (1). 

Using eq.(2) one can easily obtains the following represantation for the approximate 
solution of the initial Schrodinger equation 

H\*)=E\V), (3) 
H = h + V + H A (4) 
\*)=e g(E)\ Xo )\B). (5) 

Here Q{E) = (E — H)' 1 , H = h + Ef=i^. The state \B) = (xol*), m ore explicitly 
written as (R(xo\^) is a new unknown function, which by definition depends only on the 
relative distance R between the meson and the center of mass of the nucleonic subsystem. 



From the structure of the hamiltonian H it is clear, that the Green function G(E) 
depends on R, as a dynamical variable, and parametrically on the internuclear distances. 
This observation means, that expression (5) for the total solution never can be presented 
by the widely used RGM ansatz 

|*«?> = IXo>|S>, (6) 
Only with additional approximation of type 

(Xo\Q(Y)\ X o) « (Xo\G(V)\xo), (7) 

where V = (xol^lxo), one can reduce the general presentation (5) to the particular RGM 
form (6). 

Projecting eq.(5) onto \xo) one can immediately come to the one-dimensional integral 
equation for the unknown function \B). However, the kernel of this sort of equation can be 
calculated analytically only for the separable type of meson- nucleon potentials V^. Since 
we work with local potentials of Yukawa type, it is more practical to apply a different 
procedure. For this purpose, let us introduce the set of \xi) 1 auxiliary orthonormal square 
integrable functions defined in nuclear subspace. As the first term of this set let us take 
the wave function of the ground state of nuclei A He, \xo) 

H A \Xo)=e \xo)- (8) 
Inverting Green function in eq.(5) and using the expansion 

N 

l*> = £IAHx*>, (9) 

one can obtain the system of one-dimensional differential equations for the functions \Bi) 

N 

^{Xkln-Elx^Bi) = -e \B k )5 k0 . (10) 

i=0 

It is clear that if one restricts oneself only to the first equation of the system (10), the 
folding model for interaction of meson with nucleus will appear. 

Now let us present the results of calculation of binding energy of the ?7 c -meson with 
the nucleus A He. As it was suggested in [1], meson-nucleon interaction can be described 
by the attractive Yukawa potential 

e~ Mr 

V mN (r) = -V 

r 

with parameters varying in the following intervals: 

0.4 <V < 0.6 
0.4 < fi < 0.Q(GeV) (12) 

In accordance with the results of [2], the lowest approximation, corresponding to the 
folding model (solid line on Fig.[l]), produces the bound state in the system. It is clear 

1 For details see the Appendix 



(11) 



from the Fig. [1] that there is rather a fast convergence of the suggested procedure for all 
values for strengths of the meson-nucleon potential. In other words, it means that Green 
function 

g{E) = {E-H)-\ (13) 

which describes a relative motion of meson and the center of mass of the nucleus, can be 
reasonably well represented for the eigenvalue problem by a finite matrix of low dimension. 

The comparision of numerical values of the binding energy in the present approach and 
in the folding model [2] gives E ^^ = ™ G(N) = -5MeV, E^ e ^ =0 - &GeY) = -7.5MeV. 
Underestimation of binding in the folding model is due to neglecting in this model of 
rescattering of meson by nucleons of the target. In our formulation, this effectes fully 
taken into account in the Green function G{E). 

In [1] the binding energy we are interested in appeared to be equal to E P \ 1C11 = — 143MeV 
which is probably related with a not fully adequate A-dependence of the meson-nucleus 
potential, postulated by the authors. 



Appendix 



We will build the set of auxiliary orthonormal functions \x%) in a usual manner. First 
of all we will take a set of gaussian functions \x%) 

N-l 

l*(«0> = Il(-) 3/4 e-^ 2 , (14) 
j'=i 71 

where N is a number of the particles in the system, r« is a jacobian coordinates in the 
nuclear center of mass frame. As it is obvious this functions have a unit norm, at is a set 
of coefficients {ccji...aijjv_i}. This coefficients are chosen in the following manner, a is 
chosen so, that (xo|?" 2 |xo) — r m sri where r msr is a mean square radius of the corresponding 
nuclea A. For example, if N — 4, then this coefficients would be 

1 

"01 



2(r 2 ) 

V msr I 
1 



msr I 



"03 



(r 2 ) 

\ msr I 



Coefficients a, are chosen recursively a i+ i = «i*C, where C is some numerical coeffiecient 
chosen in such a way, that each following function \xi+i) would be wider than the previous 
one. 

Now functions \x%) are obtained in a Gramme-Schmidt recursive procedure 

IXi) = IXi) 

I v_ \Xi) ~ J2j=l{Xi\Xj) \Xj) _s 

lXl) " (i-EUfelx,) 2 ) 1 / 2 ' 1 ' 



Yukawa potential between meson and the j-th particle of the target is taken in the form 

V(Pj) = v — 
Pj 

N 

Pi = |P + H A ^l' ( 16 ) 
i=i 

where p is a distance between meson and nuclei center of mass, coefficients \j specifies 
position of the corresponding j-th particle. One could see that the following matrix 
elements takes form 

£ t M 2 



V f°° e 4(t+7 kij) 
(Xi{ a i)\ V {pk)\Xj{ a i)) = ttj= / dt 



Jo (t + 7 Wj ) 

\2 

lk .. = — . (17) 

To obtain formulas for (Xi\V\Xj) it's very convinient to use following recursive formulas 

1 j ~ 1 

(Xi\V\Xj) = (Xj\V\Xi) = -. ^j-i ,~ i — r^{(Xi\V\Xj) - ^2(Xk\Xj)(Xi\V\Xk)) 



1 l - x 

(Xi\V\Xj) = 7j _ — I- ^TE ((Xi\V\Xj) - Y,(Xk\Xi)(Xj\V\Xk))- (18) 
I 1 2^fc=i\Xfc|Xi/ ) k=i 



Due to the first equality we could always assume that i < j. 

Results for the binding energy of r) c A He system are presented on the figure below. 
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